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Measurements,  including  all  Reynolds  gtremses 
and  triple  products,  have  been  made  in  three 
turbulent  flows  with  Imbedded  longitudinal 
vortices  generated  by  skewing  of  the  mean  flow 
(Secondary  flow  of  the  first  kind*) 

(i)  flow  in  an  idealized  wing-body  junction^ 


(li)  an  Isolated  vortex  in 
boundary  layer  *=>  / 


(iii) 

layer* 


a  vortex  pair  in 


two-dimensional 


turbulent  boundary 


Selected  results  are  presented:  spatial 
variations  of  eddy  dif fusivities  and  dimensionless 
modelling  parameters  are  seen  to  be  complicated* 


1 •  Introduction 

Sufficiently  strong  lateral  deflection 
(“skewing")  of  the  mean  streamlines  in  a  shear 
layer  can  lead  to  the  generation  of  discrete 
longitudinal  vortices*  The  mechanism  is  basically 
invlscid,  given  the  presence  of  an  initial  shear: 
this  "skew- induced"  vorticity  should  be 
distinguished  from  the  generation  of  streamwlse 
vorticity  by  Reynolds  stresses,  which  is  important 
only  in  very  long,  straight  streamwlse  corners  and 
is  Confined  In  practice  to  non-circular  ducts*  In 
most  cases  of  practical  Interest  the  boundary 
layer  and  the  imbedded  vortex  will  be  turbulent 
but  the  decay  of  the  vortex  under  the  action  of 
Reynolds  stresses  is  slow,  because  its  circulation 
is  reduced  only  by  the  effect  of  the  spanwlse 
component  of  skin  friction.  That  is,  the  skew 
induced  vortices  can  Influence  the  flow  for  a  very 
long  distance  downstream.  The  present  paper  is  a 
description  of  measurements  of  the  mean  and 
turbulent  properties  of  decaying  vortices  and 
vortex  pairs:  the  process  of  formation  (being 
"invlscid")  does  not  require  such  detailed  study 
as  the  slow  decay*  The  configurations  used  are 
idealizations  of  those  found  In  practice,  but  the 
data  should  be  useful  for  developing  and  testing 
calculation  methods  Intended  for  real  life  cases* 

The  most  obvious  case  in  which  lateral 
deflection  of  a  shear  layer  leads  to  the 
generation  of  discrete  vortices  Is  divergence  of 
the  boundary  layer  flowing  over  a  surface  around  a 
tall  obstacle  protruding  from  the  surface*  The 
horseshoe  vortex  wrapped  round  the  junction 
between  a  wing  and  a  body,  an  axial  flow 
turbomachine  blade  and  a  hub,  or  a  hull  and  a  fin, 
is  formed  In  this  way*  A  single  vortex  trails 
downstream  in  each  wing-body  corner,  while 
downstream  of  the  trailing  edge  the  two  vortices 
merge  to  form  a  pair  with  the  "common  flow" 


— between  them  directed  towards  the  *  surface. 
Vortices  can  also  be  generated  on  nearly-flat 
surfaces,  by  sufficiently  strong  crossflow*  The 
vortex  pair  on  the  leeward  side  of  a  body  of 
revolution  at  incidence,  and  the  bilge  vortices 
generated  near  the  bow  of  a  ship  with  a  full  hull 
form,  are  two  examples*  A  vortex  pair  can  be 
formed  in  a  sufficiently  strong  convergent 
crossflow,  as  in  the  bottom  of  a  "S-bend"  dorsal 
Intake  or  the  wall  of  a  wind  tunnel  contraction: 
in  this  case  the  "common  flow"  is  away  from  the 
surface  and  leads  to  rapid  thickening  of  the 
boundary  layer*  The  Taylor-Cortler  vortices  that 
form  in  laminar  or  turbulent  boundary  layer  flow 
over  concave  surfaces  have  some  points  of 
similarity  with  the  decaying  vortices  discussed 
above,  but  will  not  be  discussed  further  In  the 
present  paper:  detailed  investigations  are  being 
carried  out  by  the  present  authors #  colleagues* 


The  experiments  to 
present  paper  relate  to: 


described  in 


(1)  one  leg  of  a  horseshoe  vortex  in  a  90  deg 
streamwlse  corner 

(11)  an  Isolated  vortex  in  the  boundary  layer 
on  a  flat  surface,  the  vortex  being  artificially 
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Figure  1.  Test  rig  (blower  not  shown). 
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Results 


(111)  vortex  pairs  on  flat  surfaces:  again,  the 
vortices  are  artificially  generated  far  upstream 
and  the  "common  flow'*  can  have  either  sign* 

2*  Test  rigs* 

All  the  experiments  were  done  in  a  30  in  x  5 
in  (762  mi  x  127  mm)  open  circuit  blower  wind 
tunnel  (Flg*l).  Working  section  lengths  up  to 
about  3  m  could  be  used*  The  non  Inal  tunnel  speed 
for  all  experiments  was  100  ft/sec  (30  m/s). 
Standard  Pitot  tubes,  three-hole  Conrad  yawmeters, 
and  conventional  cross  hot  wire  probes  were  used 
for  the  measurements*  The  fluctuating  signals 
from  hot  wires  were  recorded  on  analogue  magnetic 
tape,  and  later  transcribed  to  digital  magnetic 
tape  for  computer  analysis,  including 
linearization*  Statistics  involving  both  v  and 
w  components  were  deduced  from  measurements  with 
the  cross  wire-probe  rotated  through  various 
angles  about  Its  axis* 

The  corner  vortex  was  generated  by  an 
idealized  wing,  spanning  the  5  in  height  of  the 
tunnel  with  its  leading  edge  57  In*  from  the 
working  section  entrance  (Fig.2).  The  model  had  a 
ha  If -elliptical  nose  6  in.  long,  followed  by  a  slab 
of  2  in  constant  thickness*  Thus,  pressure 
gradients  were  negligible  except  in  the  leading 
edge  region,  and  there  was  no  tendency  for  the 
vortex  to  drift  away  from  the  corner  under  the 
influence  of  pressure  gradient.  The  boundary 
layer  on  the  "body"  (actually  the  floor  of  the 
wind  tunnel)  was  about  1  in. thick  at  the  leading 
edge. 


The  isolated  vortices  and  vortex  pairs  on 
flat  surfaces  are  generated  by  half-delta  wings 
mounted  in  the  settling  chamber  of  the  wind 
tunnel,  as  shown  In  Fig*  1.  As  the  flow  passes 
through  the  contraction,  the  circulation  around 
the  vortex  is  conserved,  but  the  percentage 
velocity  decrement  in  the  wake  of  the  delta-wing 
is  very  much  reduced,  so  that  at  exit  from  the  9 
to  1  two-dimensional  contraction  we  have  a 
concentrated  vortex  with  nearly  uniform  axial 
velocity*  (Although  the  contraction  is 
two-dimensional,  the  vortex  rapidly  recovers  its 
circular  shape.)  For  the  vortex  pair 
experiments,  two  delta-wing  vortex  generators,  set 
close  together,  are  used:  the  vortex  generator 
configurations  were  developed  by  flow 
visualization  experiments  in  another  wind  tunnel 
before  quantitative  measurements  began. 
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Plgur*  2.  Wing-body  junction  t*.t  section. 


Some  features  are  common  to  the  results  of 
all  the  experiments*  The  circulation  around  the 
Isolated  vortices,  and  even  the  circulation  around 
one  vortex  of  a  pair,  decreases  only  very  slowly 
with  downstream  distance.  Even  vortices  that  are 
initially  very  small  grow  quite  rapidly  to  fill 
the  boundary  layer,  the  growth  rate  being 
comparable  with  that  of  plumes  of  passive 
contaminants.  The  diameter  of  the  corner  vortex 
remains  roughly  equal  to  the  local  thickness  of 
the  body  boundary  layer  (which  is  thicker  than  the 
wing  boundary  layer)  but  the  region  of  velocity 
defect  in  the  isolated  vortices  and  vortex  pair 
projects  well  outside  the  undisturbed  boundary 
layer.  Another  feature,  common  to  all  the 
configurations,  is  the  complicated  behaviour  of 
the  various  components  of  eddy  viscosity,  which 
become  negative  in  significant  regions  of  the 
flow:  this  implies  that  the  use  of  a  simply 
behaved  eddy  viscosity  to  predict  the  spreading 
rate  of  the  vortex  Is  likely  to  be  unsatisfactory 
unless  only  gross  properties  are  required* 
Unfortunately,  the  behaviour  of  the  triple 
velocity  products  which  appear  In  the  Reynolds 
stress  transport  equations  is  also  complicated,  so 
that  although  calculation  methods  based  on 
transport  equations  are  needed  in  principle  for 
the  calculation  of  vortex  flows,  they  may  be 
difficult  to  develop  in  practice.  Although  this 
paper  is  not  directly  concerned  with  calculation 
methods,  it  should  be  pointed  out  that  a 
calculation  method  which  has  been  optimized  for 
predicting  secondary  flows  of  the  second  kind  (the 
"stress  Induced"  secondary  flows  found  in  long 
non-circular  ducts)  will  not  necessarily  be  very 
good  at  predicting  secondary  flows  of  the  first 
kind  (the  skew  induced  flows  discussed  In  the 
present  paper). 

3*1  The  corner  vortex 

This  work  Is  fully  described  by  Shabaka  (1) 
and  a  synoptic  Is  given  by  Shabaka  and  Bradshaw 
(2).  The  results  are  archived  on  magnetic  tape  in 
the  The rroosc lences  Division,  Stanford  University, 
as  part  of  the  data  set  for  the  1980/81 
Conferences  on  Complex  Turbulent  Flows  (3).  The 
"wing"  used  is  shown  In  Fig*  2,  which  also  gives 
details  of  the  initial  boundary  layers. 
Measurements  were  not  made  In  the  leading  edge 
region  in  which  the  horseshoe  vortex  is  generated, 
and  the  measurements  made  at  the  start  of  the 
parallel  section  of  the  wing  give  the  initial 
conditions  from  which  the  vortex  decays. 

Mean  velocity  contours  at  three  streamwlse 
positions  are  shorn  in  Fig.  3*  The  position  of 
the  vortex  centre,  deduced  from  results  presented 
below,  is  shown  approximately  as  a  guide  to 
understanding.  The  distortion  of  the  mean 
velocity  contours  by  the  vortex  is  relatively 
mild,  but,  as  shown  In  Fig.  A,  crossflow  angles 
of  several  degrees  are  found  near  the  wing  and 
body  surfaces.  Figure  A  hints,  and  Fig.  5  shows 
plainly,  that  the  vortlclty  is  not  a  maximum  at 
the  position  which  Fig.  A  indicates  to  be  the 
centre  of  the  vortex  but  is  largest  near  the 
surface  of  the  body,  where  of  course  it  consists 
mainly  of  the  secondary  shear  dtl/dy,  dV/dt  being 
comparatively  small.  Surface  crossflow  angles  are 
large  enough  to  be  visible  In  surface  oil  flow 


Figure  5.  Streamwise  vorticity  contours  at  x*614  j 
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Contours  of  turbulence  intensity,  of  which 
examples  are  given  In  Ref .2,  look  qualitatively 
like  the  axial  velocity  contours*  Figure  6  shows 
Isometric  views  of  the  profiles  of  the  Reynolds 
shear  stresses.  Far  from  the  corner,  the  profiles 
of  07  (the  body  shear  stress)  and  of  00  (the 
wing  shear  stress)  revert  to  two-dimensional 
behaviour,  but  07  In  particular  behaves  very 
oddly  near  the  corner,  In  a  way  that  could  not  be 
accounted  for  by  convection  of  the  basic  boundary 
layer  stress  pattern  by  the  secondary  flow.  The 
shear  stress  In  the  yz-planc,  VF,  is  fairly 
small  except  In  the  vicinity  of  the  vortex 
centre,  where  it  changes  sign  rapidly.  The  eddy 
viscosities  defined  as  -uV/dU/dy  and  -uw/dU/dz, 
have  negative  regions  near  the  wing,  and  near  the 
body,  respectively:  the  shear  stresses  in  the 
region  of  negative  viscosity  are  significant,  so 
that  the  transfer  of  axial  momentum  away  from  the 
corner  would  be  seriously  miscalculated  in  the 
region  near  the  solid  surfaces  if  an  eddy 
viscosity  based  on  two-dimensional  boundary  layer 
practice  were  deed. 

It  is  not  possible.  In  a  short  paper  to 
discuss  the  rather  complicated  behaviour  of  the 
rather  large  number  of  relevent  triple-products. 
Fig.  7  shows  profiles  of  uv?,  the  rate  of 
turbulent  transport  of  07  (the  body  shear  stress) 
in  the  y  direction  (away  from  the  body).  For 
large  values  of  s,  results  are  typical  of 
isolated  boundary  layers,  but  the  profiles  become 
greatly  distorted  nearer  the  corner,  and  uv? 
appears  to  asymptote  to  non-sero  values  at  y 
within  the  wing  boundary  layer,  although  values  at 
vry  large  y  (i.e.  in  the  two-dimensional  wing  ' 
boundary  layer  far  from  the  body)  must  tend  to 
sero. 


Figure  6.  Reynolds  stress  contours  at  x«614  mm 
3.2  Single  Isolated  vortex 

Figure  8  shows  the  axial  velocity  contours 
just  downstream  of  the  start  of  the  working 
section,  where  the  floor  boundary  layer  is  about  4 
mm  thick.  The  height  of  the  vortex  centre  is 
expected  to  be  at  about  ym\0tm9  and  Its 
s  coordinate  is  obviously  about  23mm  (the  s 
origin  being  effectively  arbitrary).  Figure  9 
shows  axial  velocity  contours,  and  secondary  flow 
vectors  about  1.3  i  downstream  of  the  start  of  the 
working  section.  For  larger  negative  values  of 
s  than  those  shown  here,  the  velocity  contours 
asymptote  to  the  approximately  two-dimensional 
state  shown  for  large  positive  s.  The  vortex 
centre  Is  quite  clearly  defined,  at  roughly  the 
same  s  position  as  In  Fig.  8,  but  at  about  20 
mm  from  the  surface,  indicating  that  the  vortex 
moves  away  from  the  surface  as  the  boundary  layer 
thickens.  (Contours  at  x»1.94  m  are  nearly 
geometrically  similar  to  those  shown  here.) 


Figure  7.  Profiles  of  triple  product  uv2  at  x*614  . 
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Figure  6.  Mean  velocity  contours  in  single 
isolated  vortex  at  x*112  wm 


As  in  the  case  of  a  corner  vortex,  maximum 
crossflow  angles  and  maximum  streamwlse  vortlclty 
occur  near  the  surface*  However,  the  region  of 
total  pressure  deficit  in  the  vortex  extends  out 
to  about  twice  the  thickness  of  the  undisturbed 
boundary  layer :  since  the  initial  total  pressure 
deficit  in  the  vortex  (Fig*  0)  is  quite  snail, 
this  low  total  pressure  fluid  oust  have  been 
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Figure  9*  Mean  velocity  contours  and  secondary 
flow  vectors  at  x“1332  n. 

transported  from  near  the  surface  by  the  secondary 
flow,  assisted  by  turbulent  diffusion* 

The  effect  of  x-vlee  convection  of  fluid  by 
the  vortex,  In  an  action  broadly  similar  to  the 
rolling  of  sheet  metal,  is  simply  demonstrated  by 
the  behaviour  of  the  skin  friction  coefficient 
shown  in  Fig.  10:  velocity  profiles  near  the 
surface  obey  the  logarithmic  law  at  all  values  of 
z,  so  that  the  skin  friction  is  roughly 
proportional  to  the  square  of  the  velocity  at, 
say,  y*  0.14.  The  skin  friction  profiles  are 
qualitatively  exploitable  by  the  transfer  of  high 
momentum  fluid  towards  the  surface  In  the  region 
of  z»-40,  and  transfer  away  from  the  surface  In 
the  region  of  z-0  :  the  secondary  peak  in  cf  at 
z*20,  corresponding  to  the  dip  In  the  velocity 
contours  near  this  position,  is  more  difficult  to 
explain*  The  secondary  velocity  contours  In 
station  9  show  some  hints  of  s  second  vortex,  with 
its  axis  at  about  z»20,  but  it  Is  not  certain 
exactly  how  this  second  vortex  Is  produced* 
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Figure  10*  Skin  friction  in  singlo  Isolated 
vortex* 
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Figure  11*  Longitudinal  component  intensity 
contours* 


Figure  12.  Primary  shear  stress  contours* 


Figure  11  shows  contours  of  longitudinal 
turbulence  intensity*  As  in  the  case  of  the 
corner  vortex,  the  intensity  contours  look  roughly 
similar  to  the  velocity  contours »  and  this 
similarity  extends  to  the  v  and  w  component 
Intensities  as  well* 

The  37  contours  shown  in  Fig.  12  are  chosen 
so  that  the  outermost  approximately  coincides  with 
the  edge  of  the  undisturbed  boundary  layer,  but  it 
Is  clear  that  the  region  of  significant  37 
within  the  vortex  is  very  much  smaller  than  that 
outlined  by  the  outermost  velocity  contour, 
bodily  rotation  of  the  fluid  In  the  vortex,  so 
that  eddies  with  a  contribution  to  37  are 

rotated  so  as  to  contribute  to  (37  Instead,  must 
be  responsible.  Figure  13  shows  that  39  (zero 
in  a  two-dimensional  flow)  is  Indeed  complicated. 
The  very  rapid  s-wlse  gradient  in  the  region  of 
*»i0  mm  lines  up  with  the  position  of  the  "tongue1* 
in  uv  shown  in  the  previous  Fig*  Contours  of 
V'  snd  "0-"  srs  shown  sepsrately  for  clarity;  in 
fact  the  zero  contour  has  several  cusps*  The 
region  of  large  positive  39  at  t«*25,  where 
u7  fails  rather  rapidly  to  zero.  Is  In  the  right 
sense  to  be  explained  as  a  rotation,  by  the  main 
vortex,  of  eddies  originally  contributing 

negatively  to  3V.  The  concentrated  region  of 
negative  OF  near  t-15  is  more  difficult  to 

explain,  as  is  that  near  *--85.  Figure  14  shows 
79,  again  with  H0fM  snd  "0-"  contours  Indicated. 
The  general  shape  of  the  contours  la  believed  to 
be  trustworthy,  but  the  large  negative  peak  at 
y«8,  i*15  appears  Implausible*  Eddy  viscosities 

corresponding  to  these  Reynolds  stresses  have  not 
yet  been  evaluated  quantitatively,  but  even  the 
casual  observation  that  07  ia  avarywhare  of  the 
same  sign,  whareaa  tha  velocity  gradient  dll/dy 
raveraas,  clearly  indies taa  that  the  eddy 
viscosity  will  be  ill-behaved. 

Figure  13  shows  profiles  of  the  triple 
product  uv7*  As  In  Fig*  7,  large  departures 
from  the  two-dimensional  behaviour  exemplified  by 
values  at  large  a  are  seen  to  occur* 
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Figure  13.  Secondary  .hear  acre.,  contour. 
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Figure  I*.  Croat  .actional  ahaar  .tr...  contour. 


Figure  15.  Profiles  of  triple  product  uv* 
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Analysis  of  the  results  of  this  experiment  is 
still  in  progress.  It  extends  the  mean- flow 
measurements  of  Tanaka  and  Suzuki  (4)  and  the  old 
NPL  results  quoted  by  Pearcey  (see  Ref .5). 

3.3  Vortex  pair  on  flat  surface  ("common  flow" 
outwards) 

Figure  16  and  17  show  mean  velocity  contours 
for  this  case.  Figure  16  shows  measurements  near 
the  front  of  the  working  section*  the  large 
scatter  In  the  points  for  the  U/Ue-0.995  contour 
being  mainly  a  consequence  of  the  complicated 
total  pressure  pattern  In  the  residue  of  the  wake 
of  the  twin-delta  vortex  generators  In  the 

settling  chamber.  In  Fig.  17*  the  outermost 
contour  Is  smooth*  but  the  contour  for  U/Uc-0.975 
is  irregular*  with  some  suggestion  that  the 
Irregularities  are  real  rather  than  just  scatter 
In  a  region  of  small  gradients.  In  general*  the 
contours  In  these  figures  are  plausible  and 

accurately  symmetrical*  albeit  about  a  value  of 
z  slightly  different  from  zero.  Figure  18  shows 
skin  friction  values  at  different  streamwise 

stations.  Again*  the  curves  are  accurately 
symmetrical*  the  dip  at  the  centreline  being  the 
obvious  result  of  Inflow  of  low  momentum  fluid 
from  positive  and  negative  z  •  The  double 

minimum  at  the  smallest  value  of  x  * 
corresponding  to  the  dip  In  the  velocity  contour 
nearest  the  wall  (Fig.  16)  Is  evidently  a  result 
of  the  disturbances  caused  by  the  vortex 
generators  because  only  a  single  minimum  Is  found 
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Conclusions 


Even  In  their  present  incomplete  state,  the 
results  of  these  experiments  on  longitudinal 
vortices  Imbedded  In  turbulent  boundary  layers 
call  attention  to  many  Interesting  features  of  the 
flows.  In  which  the  qualitative  obviousness  of  the 
mean  velocity  contours  conceals  a  great  deal  of 
non-obvlous  turbulence  behaviour,  which  will  have 
to  be  better  understood  If  the  mean  velocity 
contours  are  to  be  predicted  quantitatively  • 
Longitudinal  imbedded  vortices  are  such  a  common 
feature  of  three-dimensional  turbulent  flows  that 
serious  attempts  to  develop  a  calculation  method 
are  long  overdue*  The  results  indicate  that 
methods  using  eddy  viscosity  or  eddy  dlffusivlty 
are  unlikely  to  be  satisfactory,  but  even 
transport  equation  methods  cannot  be  expected  to 
reproduce  all  the  details  of  the  complicated 
profiles  of  Reynolds  stresses  shown  above.  The 
simple  configurations  tested  In  the  present 
experiments  correspond  to  a  rather  wide  variety  of 
real  life  flows,  of  importance  in  such  diverse 
fields  as  turbomachinery,  aeronautical 
aerodynamics,  ship  hydrodynamics,  and  Internal 
flows. 
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